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SUMMARY 

A general method of design is developed for two-dimensional 
unbranched channels with prescribed velocities as a function of arc 
length along the channel walls. The method is developed for both 
incompressible and compressible; irrotational, nonviscous flow. Two 
types of compressible flow are considered: the general type, with the 

ratio of specific heats y equal to 1.4, for example; and the linear- 
ized type m which x 1S equal to -1.0. The design method gives com- 
plete information concerning the flow throughout the channel. 

Five numerical examples are given including three elbow designs 
with the same prescribed velocity as a function of arc length along the 
channel walls but with incompressible, linearized compressible, and 
compressible flow. It is concluded that if a nonviscous gas with 
arbitrary y (1.4, for example) were to flow through a channel designed 
for linearized compressible flow (y = -1.0) > the resulting velocity 
distribution along the channel walls would be nearly the velocity dis- 
tribution prescribed for the linearized compressible flow, at least if 
the linearized flow were selected so that the densities are equal for 
both types of flow at the maximum and minimum velocities and if the 
ratio of these velocities is not too large (2:1 m the numerical 
examples) . 


INTRODUCTION 

There are two general types of theoretical problem m two- 
dimensional fluid motion: (l) the direct problem, in which the distri- 

bution of velocity is determined for a prescribed shape of boundary, 
and (2) the inverse problem, in which the shape of boundary is determined 
for a prescribed distribution of velocity along the boundary. The direct 
problem is an analysis problem; the inverse problem is a design problem. 
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This report is concerned with the inverse, or design, problem for two- 
dimensional, irrotational flow in unbranched channels with prescribed 
velocities as a function of arc length along the channel walls. 

The design of channels with prescribed velocities is important 
because: (l) Boundary- layer separation losses can be avoided by pre- 

scribed velocities that do not decelerate rapidly enough to cause 
separation, (2) shock losses in compressible flow and cavitation 
in incompressible flow can be avoided by prescribed velocities that do 
not exceed certain maximum values dictated by these phenomena, and 
(3) for compressible flow the desired flow rate can be assured by pre- 
scribed velocities that do not result in M choke flow" conditions. 

Several methods of channel design have been developed for particular 
application (references 1 and 2, for example) . In reference 1 a design 
method is developed for accelerating elbows m which the velocity 
increases monotonically along the channel walls. The method is developed 
for incompressible and linearized (y = -1.0) compressible flow. The 
velocity distribution along the channel walls is not arbitrary and the 
design method applies to elbows only. In reference 2 a design method is 
developed for straight, symmetrical channels with contracting or expand- 
ing walls. The method is developed for incompressible flow and the 
velocities are prescribed not as a function of arc length along the 
channel walls but as a function of circle angle in the transformed circle 
plane. A more general design is suggested in reference 3 but no attempt 
is made to develop and apply the method. 

In the present report a general method of design is developed for 
two-dimensional, unbranched channels with prescribed velocities as a 
function of arc length along the channel walls. The method is developed 
for both compressible and incompressible, irrotational, nonviscous flow 
and applies to the design of elbows, diffusers, nozzles, and so forth. 

Two types of compressible flow are considered: the general type with 

arbitrary value of y (1.4, for example) and the linearized type with 
y equal to -1.0. In general, if the prescribed velocity along one 
channel wall differs from that along the other, the channel turns so 
that the downstream flow direction is different from the upstream 
direction. This change in flow direction cannot be arbitrarily chosen 
but depends on the prescribed velocity distribution along the walls. 
Equations are developed for computing this change in flow direction for 
an arbitrary prescribed velocity distribution with incompressible or 
linearized congress lble flow. Two methods of solution have been devel- 
oped for the design method and are presented m separate reports. In 
this report (part I) solutions are obtained by relaxation methods 
(reference 4) . This method of solution results m complete information 
concerning the distribution of flow conditions throughout the channel 
and, in addition, can be used to obtain nonlinear solutions for com- 
pressible flow with arbitrary values of y. In reference 5 (part II) 
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solutions are obtained by means of a Green* s function. This method of 
solution is limited to incompressible and linearized (y = -1.0) com- 
pressible flow, but the method is more rapid than relaxation methods, 
provided information within the channel is not required. 

The design method reported herein was developed at the MCA Lewis 
laboratory during 1950 and is part of a doctoral thesis conducted with 
the advice of Professor Ascher H. Shapiro of the ^ssachusetts Institute 
of Technology. 


THEORY OF DESIGN METHOD 

The design method is developed for two-dimensional channels with 
prescribed velocities along the channel walls. The prescribed velocity 
is arbitrary except that stagnation points (zero velocity) cannot be 
prescribed. This exception limits the design method to unbranched 
channels . 


Preliminary Considerations 

Assumptions. - The fluid is assumed to be nonviscous and either 
compressible or incompressible. The flow is assumed to be two- 
dimensional and irrotational. 

The assumption of two-dimensional, nonviscous, irrotational motion 
limits the design method m practice to channels with thin (negligible) 
boundary layers, such as exist near the entrance to the channel or after 
a rapid acceleration of the flow through a contraction m the channel. 
Even if the boundary layer is thin, the design method is limited to (and 
finds its most useful application for) prescribed velocity distributions 
that, from boundary- layer theory, do not decelerate fast enough to 
result m separation of the boundary layer, which separation alters the 
"effective" shape of the channel and completely changes the character of 
the flow. 

In some channels with fully developed turbulent boundary layers the 
design method might be expected to yield results that are satisfactory 
(although approximate) because for this type of flow the rotational 
motion occurs primarily m the regions close to the channel walls. In 
channel walls with thick or fully developed laminar boundary layers the 
design method cannot be used, because not only is the rotation of the 
flow important m most of the channel but, if the channel bends, 
important secondary flows develop that are not considered by the two- 
dimensional design method. 
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Flow field. - The flow field of the two-dimensional channel is con- 
side redTTcTTii - in the physical xy-plane where x and y are Cartesian 
coordinates expressed as ratios of a characteristic length equal to the 
constant channel width downstream at infinity. (All symbols are defined 
in appendix A.) 

At each point m the channel (fig. 1) the velocity vector has a 
magnitude' Q and a direction 9 where Q is the fluid velocity 
expressed as the ratio of a characteristic velocity equal to the con- 
stant channel velocity downstream at infinity. For convenience, the 
velocity Q is related to a velocity q by 

q = Qla (!) 

where q is the velocity expressed as a ratio of the stagnation speed 
of sound and the subscript d refers to conditions downstream at 
infinity. 

The flow direction 6 at each point m the channel is measured 
counterclockwise from the positive x-axis. From figure 1 

dx = ds cos 9 (2a) 

dy = ds sm 9 (2b) 

where ds is a differential 'distance in the direction of Q, that is, 
along a streamline . 

Stream function and velocity potential . - if the condition of con- 
tmuity is satisfied a stream function 'p can be defined such that 

d\p = P Q dn (3) 

where p is the fluid density expressed as the ratio of a characteristic 
density equal to the stagnation density and where dn is a differential 
distance measured normal to the direction of Q, that is, normal to a 
streamline. Along a streamline, dn is zero so that from equation (3) 
the stream function \[/ is constant. 

If the condition of irrotational fluid motion is satisfied a 
velocity potential cp can he defined such that 

dcp = Q ds (4) 

Normal to a streamline, ds is zero so that from equation (4) the 
velocity potential cp is constant. Thus lines of constant cp and ^ 
are orthogonal m the physical xy-plane. 
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Outline of method. - Solutions for two-dimensional flow depend on 
known conditions imposed along the boundaries of the problem. In the 
inverse problem of channel design the geometry of the channel walls m 
the physical xy-plane is unknown. This unknown geometry apparently 
precludes the possibility of solving the problem m the physical plane 
and necessitates the use of some new set of coordinates , that is, a 
transformed plane, in which to solve the problem. These new coordinates 
must be such that the geometric boundaries along which the velocities 
are prescribed are known in the transformed plane. It is also desirable 
for mathematical simplicity, that the coordinate system in the trans- 
formed plane be orthogonal m the physical plane. A set of coordinates 
that satisfies these requirements is provided by cp and \|/ , which are 
orthogonal m the physical xy-plane and for which the geometric bound- 
aries are known constant values of ^ m the transformed cp\J/-plane. 

The distribution of velocity as a function of cp along these boundaries 
of constant \j / is known because, if 

Q = Q(s) 

is prescribed, equation (!) integrates to give 

cp = cp(s) 

From which equations, 

Q - Q(cp) 

The technique of the proposed method of channel design is therefore 
to obtain a differential equation for the distribution of velocity in 
the cp\|/-plane. The velocity distribution obtained from the solution of 
this equation is then used to obtain the distribution of flow direction, 
from which distribution the channel walls m the physical xy-plane are 
obtained directly. The differential equation for the distribution of 
velocity m the Cpty -plane is nonlinear (for compressible flow with j 
other than -1.0) and is solved by numerical methods (relaxation methods) 


Differential Equation for Distribution of Velocity 
in Transformed cp\J/- Plane 

The differential equation for the distribution of velocity in the 
transformed cp^ -plane is obtained from the equations for continuity and 
irrotational fluid motion expressed in terms of the transformed coordi- 
nates cp and ^ . 
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Continuity . - The continuity equation expressed in terms of cp and 
becomes (appendix B) : 


l/5logep & log e Q \ 
p ^ + 5cp J + "Sty 


0 


(5) 


Irrotational fluid motion. - The equation for irrotational fluid 
motion, expressed m terms of cp and ijr, becomes (appendix B) : 


5 lp ge ft _ ae = 

S\j r 3cp 


( 6 ) 


Differential equation for distribution of velocity . - The second 
order partial differential equation for the distribution of log e Q 
in the transformed cpnjr-plane 1S obtained by differentiating equations (5) 
and (6) with respect to cp and i(r, respectively, and combining to 

eliminate Thus, 

5 2 log e p a 2 log e Q d log e p fb log e p 5 log e 

+ i? \ ^P~ + / 

2 

2 d log e Q d log e p a log e Q 

p —3? 3*— + o — rr- - 0 < 7 > 

Equation (7) , together with a relation between p, Q, and q^, determines 

the distribution of log e Q m the cp^Jr-plane for compressible flow with 
a given value of q^ and for arbitrarily prescribed variations in 
log e Q along the boundaries of constant \|r. 

Density . - The density p is related to the velocity q by 
(reference 6, P« 26, for example) 

p = (l - ^ I 2 ) r_1 (8a) 

which, from equation (l), becomes 

1 

p = f 1 " q2 q a. 2 ) ( 8b ) 
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Equation (8b) relates the density p to the velocity Q for a given 
value of q^. 


Incompressible flow . - For incompressible flow p is constant and 
equal to 1.0 so that equation (7) becomes 


d 2 log e Q cl 2 log e Q 

_ — + — = 0 


Sep* 




( 9 ) 


Equation (9) determines the distribution of log e Q, m the cp^-plane 
for incompressible flow. 


Channel Wall Geometry 

After equation (7) or (9) has been solved to obtain the distribution 
of log e Q m the transformed ep^-plane (for arbitrary variations in 
log e Q with cp along the boundaries of constant i|r), the geometry of 
the channel walls m the physical xy-plane can be determined from the 
resulting distribution of flow direction 0. 

Flow direction 8 . - The distribution of flow direction 0 along 
a streamline ( constant \|r) is obtained from equation (6), which inte- 
grates to give 

r ^ iog e q 

9 = P 5* dCp dOa) 

where the subscript ijr indicates that the integration is taken along a 
line of constant \|r and where the constant of integration is selected 
to give a known value of 0 at one value of CD along each streamline. 
The integrand m equation (10a) is obtained from the distribution of 
l°g e Q; which is known from the solution of equation (7) or (9). 

The distribution of flow direction 0 along a velocity-potential 
line (constant cp) is obtained from equation (5), which integrates to 
give 


0 = 



1 /S log e P d log e q\ 

p ( — w~ + — 


(10b) 
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where the subscript cp indicates that the integration is taken along 
a line of constant cp and where the constant of integration is selected 
to give a known value of 6 at one value of i|r along each velocity- 
potential line. As for equation (10a), the integrand in equation (lOtf) 
is known from the distribution of log e Q obtained from the solution 
of equation (7) or (9). 

Channel wall coordinates . - The variation m x along a line of 
constant \|r in the cp^-plane is given by 


dx / dx ds \ 

Sep yds dep 

which, combined with equations (2a) and (4), integrates to give 



Likewise, 



(lla) 


(lib) 


(11c) 


y 



di|r 


(lid) 


where the constants of integration are selected to give known values of 
x or y at one value of cp along each streamline or at one value of 
\|r along each velocity-potential line. Equations (lla) to (lid) deter- 
mine the distribution of x and y m the transformed cp^-plane or, 
which is the same thing, the shape of the streamlines and velocity- 
potential lines m the physical xy-plane. In particular, equations (lla) 
and (lie) when integrated along the boundaries of constant \|r m the 
cp^r-plane determine the shape of the channel walls . 


Turning angle . - In general, if the prescribed velocity distribution 
along one channel wall differs from the distribution along the other 
wall, the channel deflects an amount AS, which is the difference m 
flow direction far downstream and far upstream of the region in which 
the prescribed velocity distribution varies. In reference 5 it is shown 
that for incompressible flow the turning angle A 9 is given by 
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A9 = e d - e u 



d log e 

5^ Jl.O 




( 12 ) 


where the subscript u refers to conditions upstream at infinity and 
where the subscripts 0 and 1.0 refer to the channel boundaries along 
which i ) r equals 0 and 1.0, respectively. A similar equation will be 
given later for the case of linearized compressible flow. 


Linearized Compressible Flow 

The nonlinear differential equation (7) for the distribution of 
velocity in the cpijr-plane with compressible flow becomes linear and is 
considerably simplified if a linear variation in pressure with specific 
volume (l/p) is assumed. This linear relation between pressure and 
specific volume was first suggested by Chaplygin (reference 7) in order 
to linearize the differential equations for two-dimensional compressible 
flow m the holograph plane. 

Density . - If a linear variation m pressure with specific volume 
is assumed, the density p* is related to the velocity q* by 
(appendix C) 


where 


and 


♦ ,- 2 )- 1/2 

(13) 

= kj_p 

(13a) 

= k 2 q 

(13b) 


where the constants k^ and kg have been determined so that values 
of p given by equation (13) equal the values of p given by equa- 
tion (8a) for any two selected values of q (designated by q a and 
q b ) . Thus, 
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where p a and p^ are determined by equation (8a) for the selected 
values of q a and q-^, respectively. A discussion of the selection of 
q a and q^ is given later. It will be noted that, if x ls equal to 
-1.0, equation (8a) has the same form as equation (13). 

Stream function and velocity potential. - For the case of linearized 


compressible flow it is 
the velocity potential 

convenient to define the stream function 
cp* by . 

\|r* and 


di|r* = p*q* dn 

(15) 

and 

dCp* = q* ds 

(16) 


Continuity . - The continuity equation expressed m terms of cp* 
and becomes (appendix D) 


where 


5 l°g e u 

di|r* 


(17) 


u = 


+ \f: 


1 + q* 


(18) 
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or, conversely 


q 


* 



(19) 


Irrotational fluid motion . - The equation for irrotational fluid 
motion, expressed in terms of cp* and i|r* becomes (appendix D) 


d log e u 


be 

bw* 


0 


( 20 ) 


Differential equation for distribution of log e u. - The partial 
differential equation for the distribution of log e u in the cp *4r*-plane 
is obtained by differentiating equations (17) and (20) with respect to 


cp* and i|r*, respectively, and combining to eliminate 


a 2 9 

Sep* di|r*’ 


Thus 


S 2 log e u S 2 log e u 
S^P * 2 + S ^* 2 


( 21 ) 


Equation (21) determines the distribution of log e u m cp*^ *-plane for 
linearized compressible flow with a given value of q<p and for arbi- 
trarily prescribed variations m log e Q, related to log e u by equa- 
tions (1), (13b), and (18), along the boundaries of constant 
Equation (21) is linear and is, like equation (9) for the case of incom- 
pressible flow, the equation of Laplace. Thus an incompressible flow 
solution for the distribution of log e Q m the cp^r-plane is also a 
linearized compressible flow solution for the distribution of log e u 

m the cp*^*-pl&ne. The transformation from the cp^-plane is different, 
however, from the transformation from the cp*\|/* -plane so that different 
channel shapes result m the xy-plane. 


Flow direction 6 . - The distribution of flow direction 0 along 
a streamline ( constant \| *) is obtained from equation (20) , which inte- 
grates to give 


0 = 



S log u 
°e 


dep* 


(22a) 
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Likewise, the distribution of flow direction Q along a velocity- 
potential line (constant cp*) is obtained from equation (17), which 
integrates to give 

P d log u 

e --J ,-s?- 4 ** (2: 

'-'cp* 

Equations (22a) and (22b) for linearized compressible flow correspond 
to, and are used in the same manner as, equations (10a) and (10b) for 
the usual type of compressible or incompressible flow. 

Channel wall coordinates . - The variation m x along a line of 
constant i|r* in the cp*\|r ¥ -plane is given by 


dx _ / dx ds \ 

W = \ ds d 

which combined with equations (2a) and (16) integrates to give 


Likewise; 


x 



cos 0 

q* 


d<»* 


X 



sm 9 

p*q* 


d\|r* 


(23a) 


(23b) 


y 



sin 9 

q* 


d®* 


(23c) 


y 



cos 9 

p*q* 


di|r* 


(23d) 


Equations (23a) to (23d) determine the distribution of x and y in the 
transformed cp*i|r*-plane or, which is the same thing, the shape of the 
streamline and velocity- potential lines in the physical xy-plane. In 
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particular, equations (23a) and (23c), when integrated along the bound- 
aries of constant \Jf* m the cp*i|r* -plane, determine the shape of the 
channel walls. Equations (23a) to (23d) for linearized compressible 
flow correspond to, and are used m the same manner as, equations (lla) 
to (lid) for the usual type of compressible or incompressible flow. 

Turning angle . - In reference 5 it is shown that for linearized 
compressible flow the turning angle, or difference m flow direction 
far downstream and far upstream of the region m which the prescribed 
velocity distribution vanes along the channel walls, is given by 



d log e u 
S'?* 




u 



(24) 


where Ai(r * is the value of \|r* along the left boundary (channel wall) 
when faced in the direction of flow if the value of if* along the right 
boundary is zero and where the subscript A if* refers to the boundary 
along which \|r* is equal to 


NUMERICAL PROCEDURE 

The channel design method of this report was developed for three 
types of fluid flow: (l) compressible, (2) incompressible, and 

(3) linearized compressible. Although the numerical procedures of the 
design method are similar for each type of fluid, the procedures differ 
in detail and are therefore considered separately xn this section. 


Compressible Flow 

The numerical procedure for channel design with compressible flow 
(y = 1.4, for example) is as follows: 

(l) The velocity is specified as a function of arc length along 
that portion of the channel walls over which the velocity varies 


q = q(s) 


or q^ is specified, and 


Q = Q(s) (25) 

where s is arbitrarily equal to zero at that point along one channel 
wall where the velocity first begins to vary. 
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(2) The channel wall boundaries of the flow field m the transformed 
cp\|/-plane are straight and parallel lines of constant \|r extending 
indefinitely far upstream and downstream between CD equals where 

cp is arbitrarily equal to zero at that point on the channel wall at 
which s is equal to zero. The value of i|r along the right channel 
wall when faced in the direction of flow (direction of positive cp) is 
arbitrarily set equal to zero m which case the value of \|r along the 
left channel wall (A\(r) is obtained by integrating equation (3) across 
the channel at a position far downstream where flow conditions are 
uniform 


— P,^ (26) 

(3) The distribution of log e Q as a function of CD along the 
boundaries m the cpilr-plane is obtained by integrating equation (4) 
between limits so that 

f s 

cp = / Q ds = cp(s) (27) 

Jo 

which together with equation (25) gives the distribution of log e Q 
along the boundaries m the cp\|r-plane 

Ic^ Q = f (cp) (28) 

The integration indicated by equation (27) is carried out numerically 
for arbitrary distributions of Q as a function of s. 

(4) If the velocities prescribed along one channel wall differ 
from those along the other wall, the channel will, m general, turn 
the flow. This turning angle cannot be determined exactly for com- 
pressible flow until the channel design is completed. However, it will 
be shown that this turning angle is only slightly greater than that 
resulting for linearized compressible flow with the same prescribed 
velocity and with a suitable selection for q a and q-^ m equa- 
tions (14a) and (14b) . This latter turning angle for linearized com- 
pressible flow is given by equation (24), which can be integrated numer- 
ically for the arbitrary distribution of log e u = f (cp) corresponding 
to equation (28) . 

(5) In order to solve equation (7) for the distribution of log e Q 

in the cp^-pl a ne it is convenient to eliminate the density terms from 
equation (7) by means of equation (8b). Thus, equation (7) becomes 
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b 2 log e q b 2 log e q f b log q\ 2 / d log e q ‘ 2 

A ~^ + B + 4D ( 


= 0 


( 29 ) 


where 


A = 


l-Ciq 2 


1 T ~ 1 o 2 

1 - ~ q 


1+1 

r-i 


and 


B = 1.0 


-q 


4C = 


(i * ^ »*) 


1 - ¥ * 2 


2r_ 

r - 1 


4D = 


-q 


-] r-i 2 

1 - — q 


Although equation (29) is nonlinear, it can be solved by relaxation 
methods (references 4 and 8, for example) . A grid of equally spaced 
points, at each of vhich the value of log e Q is to be determined, is 

placed m the flow field between the channel wall boundaries. The grid 
is extended upstream and downstream sufficiently far so that constant 
values of log e Q are obtained across the channel by the relaxation 

methods. In the numerical examples to be presented six or eight grid 
spaces were used across the channel. In example III the number of grid 
spaces was reduced from eight to four with negligible effect upon the 
resulting channel design. The values of log e Q at each grid point 

were relaxed to five significant figures. If the same velocity distri- 
bution is prescribed along both walls, the channel is symmetrical so that 
the velocity distribution in only one half of the channel need be deter- 
mined by relaxation methods. 
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(6) After log e Q has been determined at each grid point m the 

cp^-plane the distribution of 0 is determined by equations (10a) and 
(10b), which are integrated numerically. The constants of integration 
m equations (10a) and (10b) are determined to give a specified value 
of 0 at one point m the channel (far upstream, for example) . The 
integrands in equations (10a) and (10b) are determined by numerical 
methods (tables I to VII, reference 4, for example) from the known 
values of p and log e Q, at each of the grid points. If it is desired 
to know the flow direction along the channel walls only, equation (10a) 
can be solved along the channel wall boundaries ^ = 0 and ^ = A^r 
only. If it is desired to know 0 everywhere in the channel, the 
recommended procedure is to determine the variation m 0 along the 
mean streamline (i|r = (Ai(r) /2) by equation (10a) and to determine the 
variation m 0 along each velocity-potential line from the previously 
determined values on the mean streamline by equation (10b) . 

(7) After the distribution of log e Q and 0 are known in the 

qpty-plane, the shapes of the streamlines and the velocity-potential lines 
m the physical xy-plane or, which is the same thing, the distributions 
of x and y m the transformed cp\|r-plane are determined by the 

numerical integration of equations (lla) to (lid) . The constants of 
integration m these equations are determined so that specified values 

of x and y occur at one point in the flow field. The recommended 

procedure is to determine the variation m x and y along the mean 
streamline by equations (lla) and (lie) and to determine the variation 
m x and y along each velocity-potential line for the previously 

determined values on the mean streamline by equations (lib) and (lid) . 

If it is desired to know the x and y coordinates from the channel 
walls only, equations' (lla) and (lie) can be solved along the channel 
wall boundaries ^ = 0 and ^ = A^r only. 


Incompressible Flow 

The numerical procedure for channel design with incompressible flow 
(p = 1) is similar to that just outlined for compressible flow, but 
with the following differences: 

(1) The velocity is specified as a function of arc length by 
equation (25) alone. The constant q^ is not considered, because it 
does not exist. 

(2) The value of ^ along the left channel wall (A\|r) is equal to 
1.0 instead of the value given by equation (26). 

(3) The distribution of log e Q as a function of cp along the 
channel wall boundaries m the cp\|r-plane is the same as that obtained 
from equations (25) and (27) and given by equation (28) . 
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(4) The turning angle A 8 of the channel is given by equation (12) . 

(5) The distribution of log e Q m the cp^Jr-plane is obtained from 
the solution of equation (9) by relaxation methods. 

(6) After log e Q has been determined at each grid point between 

the channel wall boundaries m the CD\|r-plane, the distribution of 8 is 
determined by equations (10a) and (10b) as indicated previously for 
compressible flow, but with p equal to unity. 

(7) After the distribution of log e Q and 0 are known m the 

Cp^r-plane, the shapes of the streamlines and velocity-potential lines 
in the physical xy-plane are determined by equations (lla) to (lid) as 
indicated previously for compressible flow, but with p equal to unity. 


Linearized Compressible Flow 

The numerical procedure for channel design with linearized compres- 
sible flow (y = -1.0) is similar to that previously outlined for com- 
pressible flow, but with the following differences: 

(l) The velocity q is specified as a function of arc length along 
the channel walls by q(s) or by q^ and equation (25) . For each pre- 
scribed velocity there are an infinite number of linearized compressible 

flow solutions depending on the selected values of q and q m 

a b 

equations (14a) and (14b) . However, for values of q & and q^ within 

the range of q prescribed along the channel walls (and therefore 
everywhere in the channel), the solutions, that is, channel shapes, 

probably differ only in small detail. The best solution is that most 

nearly like the nonlinear compressible solution with arbitrary value of 
X (1.4, for example). In the numerical examples of this report it is 
shown that if q a and q-^ are equal to the maximum and minimum values 

of q a good solution results, at least if the ratio of these prescribed 

velocities is not too large (2:1 m the numerical examples). On the 
other hand, if continuity is to be satisfied for a gas with the correct 
value of x (1*4:, for example) upstream and downstream of the region 
of the channel m which the prescribed velocities vary, then q a and 
q^ must equal q Q and q^. 

After q a and q-^ have been selected the velocity distribution 

q(s) is expressed as q*(s) by equation (13b) where k 2 is given by 
equation (lib) so that 
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q* = q* ( s ) (30) 

The velocity q* is then expressed as u hy equation (18) so that 

u = u(s) (31) 

In the particular case where the selected value of q a is equal to 
q-^ the value of kg is given by equation (C4b) m appendix C where 
the significance of this particular case is also discussed. 

(2) The solution is obtained m the transformed co*^* -plane where 
cp* and \J i* are defined by equations (16) and (15), respectively. If 
the value of \(r* along the right channel wall when faced m the 
direction of q* is zero, the value of ijr* along the left wall (Aijr *) 
is obtained by integrating equation (15) across the channel at a 
position far downstream where flow conditions are uniform 

At* = pj qj (32) 

(3) The distribution of log e u as a function of CD* along the 
channel wall boundaries m the cD*\(r*-plane is obtained by integrating 
equation (16) between limits similar to those discussed previously for 
compressible flow so that 


C p* = - / q* ds = cp*(s) (33) 

d) 

which together with equation (31) determines the distribution of log e u 
along the channel wall boundaries in the cp*\|r*-plane 

log e u = f(cp*) (34) 

(4) The turning angle A 6 of the channel is given by equation (24) . 

(5) The distribution of log e u in the co*^*-plane is obtained 
from the solution of equation (21) by relaxation methods. 

(6) After log e u has been determined at each grid point between 
the channel wall boundaries in the cp*^*-plane, the distribution of 0 
is determined by equations (22a) and (22b) m a manner similar to that 
outlined previously for compressible flow. 

(7) After, the distribution of log e u and 6 are known in the 
cp*^*-pl a ne, the shapes of the streamlines and the velocity- potential 
lines in the physical xy-plane are determined by equations (23a) to (23d) 
in a manner similar to that outlined previously for compressible flow. 

The velocities q* m equations (23) are obtained from the known values 
of u and the densities p* are given by equation (13) . 
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NUMERICAL EXAMPLES 

The channel design method has been applied to five examples listed 
below; 


Examples 

Type of channel 

Type of flow 

I 

Reducing section 

Incompressible 

II 

Converging section 

Incompressible 

III 

Elbow 

Incompr e s s lble 

IV 

Elbow 

Linearized compressible 

V | 

Elbow 1 

Compressible (r = 1.4) 


Example I 

The first numerical example is the design of a reducing section m 
a straight channel such that the upstream velocity is half the downstream 
velocity. The solution is for incompressible flow. 


Prescribed velocity distribution, 
function of arc length s along both 

- The prescribed velocity as 
channel walls is given by 

a 

Q = 0.5 

(s <0) 




, 2 3 

n _ 1 s s 

Q 2 + 6 ' 27 

(0 < s < 3.0) 


> 

(35) 

Q = 1.0 

(s < 3.0) 

J 




The prescribed velocity given by equation (35) is plotted m figure 2. 
Equation (35) together with equation (27) results in 


cp = 0.5 s 

s s 3 s 4 

C ° “ 2 + 18 " 108 

cp = -0.75 + s 


(s < 0) 


'N 


(0 < s < 3.0) 
(s < 3.0) 


(36) 


From equations (35) and (36) , log e Q is a known function of cp, which 
function is plotted m figure 3. 
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Results . - The results of example I are presented m figures 4 to 7. 

In figure 4 lines of constant velocity Q and flow direction 9 
are plotted m the transformed cp\[r-plane. The flow direction 9 is 
constant and equal to zero along the mean streamline (i|r = 0.5) indi- 
cating that the center line of the channel is straight. The maximum 
absolute values of 9 occur along the channel walls. The solution is 
symmetrical about the mean streamline. The lines of constant Q and 
9 are orthogonal (see appendix E) . If (&S)q is the spacing of lines 
of constant 9 measured along lines of constant Q and if (&S)g is 

the’ spacing of lines of constant Q measured along lines of constant 
9 , equation (F5) m appendix F gives 


(&S) Q 



Q = 


/2rt/l80 
\ 0.03 


277 Q 


(37) 


In figure 5 lines of constant x and y are plotted on the trans- 
formed cpty-plane. Along the mean streamline (>Jf = 0.5) the value of 
y is constant and equal to zero indicating, as before, that the center 
line of the channel is straight. The lines of constant x and y are 
orthogonal (appendix E) . The solution is symmetrical. The ratio of 
the spacing of lines of constant x and y is given by equation (F6) 
of appendix F 


_ by _ 0.2 _ 1 Q 

(SS) 6x 0.2 

so that the system of curves forms a square network. 

In figure 6 lines of constant cp and \| r (velocity potential and 

streamlines, respectively) are plotted m the physical xy-plane. * The 
shape of the channel walls is that required to result m the prescribed 
velocity distribution given by equation (35) and plotted m figure 2. 

The downstream channel width is 1.0 by definition. The upstream channel 
width is 2.0 in order that the upstream velocity be half the downstream 
velocity. As usual the streamlines and velocity potential lines are 
orthogonal (appendix E) and, for equal increments of cp and form' 
a square network (equation (F7), appendix F, with p equal to 1.0). 

In figure 7 lines of constant Q and 9 are plotted m the 

physical xy-plane. The lines of constant Q and 9 are orthogonal 
(appendix E) . The ratio of the spacing of lines of constant Q and 9 
is given by equation (F8) of appendix F 
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(SS) Q = M q 0 /2n/l80 \ a = _5_ Q 

T&Sj^ \ & Q/ \ 0.03 2.7 ^ 

which is the same as that for the same lines of constant Q and 0 in 
the cpty-plane (see equation (37)). 

Example II 

The second numerical example is the design of a converging section 
that funnels the fluid from an infinite expanse into a straight channel 
of unit width. Far upstream the channel walls are straight and converge 
at a 90° angle. The solution is for incompressible flow. 


Prescribed velocity . - The prescribed velocity as a function of 
arc length s along both channel walls is given by 


Q = 


-2 


rt(s-2) 




(s < 0) -N 

(0 < s < 4) > (38) 

(s >4) J 


Q = 1.0 

The prescribed velocity given by equation (38) is plotted m figure 8. 
Equation (38) together with equation (27) results m 


cp=^l°g e (l-f 


1 0 , 1 s 

CD = — S + 

^ J t 

qp = - 2 + s 

3n 


1/7 3\ s 3 l/2 \s 4 

2it 2 " 8 \2jt 2/ 3 + 32 Ijt / 4 


(s < 0) ^ 

(0 < s < 4)>(39) 

(s > 4) j 


From equations (38) and (39) ; log e Q is a known function of cp > which 
function is plotted m figure 9. 


Results . - The results of example II are presented in figures 10 


to 12 . 
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In figure 10 lines of constant velocity Q and flow direction 0 
are plotted in the transformed cpi|r-plane. The flow direction 0 is 
constant and equal to zero along the mean streamline (\|r = 0.5) indi- 
cating that the center line of the channel is straight. The solution 
is symmetrical about the mean streamline. As for example I the lines 
of constant Q and 0 are orthogonal. The ratio of the spacing of 
lanes of constant Q and 0 is given by equation (F5) m appendix F 



In figure 11 lines of constant cp and \|r are plotted m the 
physical xy-plane. The shape of the channel walls is that required to 
result m the prescribed velocity distribution given by equation (38) 
and plotted m figure 8. As usual the streamlines and velocity poten- 
tial lines are orthogonal (appendix E) and, for incompressible flow 
with equal increments of cp and \|r, form a square network (appendix F) . 

In figure 12 lines of constant Q and 0 are plotted m the 
physical xy-plane. The lines of constant Q and 0 are orthogonal 
(appendix E) , and the ratio of the spacing of c lines of constant Q and 
0 is the same as that given for the same lines of constant Q and 0 
in the cpi|r-plane (equation (40)). 


Example III 

The third numerical example is the design of an elbow for which 
the upstream velocity is half the downstream velocity. The prescribed 
velocities are such that no deceleration occurs anywhere along the 
channel walls. The solution is for incompressible flow. 

Prescribed velocity distribution . - Along both walls upstream of 
the elbow the velocity Q is equal to 0.5, and along both walls 
downstream of the elbow Q is equal to 1.0. The transition from Q 
equals 0.5 to 1.0 along both walls of the elbow will be the prescribed 
velocity distribution as a function of arc length given by equation (35) 
for example I and plotted m figure 2. In terms of log e Q as a 
function of cp this prescribed velocity distribution is given by equa- 
tion (36) and is plotted in figure 3. Although this velocity distri- 
bution is the same for both walls, the distribution on the outer wall 
(wall with larger radii of curvature) is shifted in the positive cp 
direction an amount equal to 2.25 relative to the distribution on the 
inner wall. Thus, a velocity difference exists on the two walls at 
equal values of cp, as shown in figure 13. The greater this difference 
in velocity and the greater the range m cp over which velocity differ- 
ences exist, the greater is the elbow turning angle. For the prescribed 
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velocity distribution given m figure 13 the elbow turning angle given 
by equation (12) was 89.37 degrees compared with a value of 89.36 degrees 
obtained from the relaxation solution. 

Results. - The results of example III are presented in figures 14 
to 16 and m tables I and II. (The numerical results for examples III, 
IV; and V axe tabulated m tables I to VI to enable a detailed comparison 
of the three elbow designs with the same prescribed velocity Q distri- 
bution as a function of arc length but with incompressible (example III); 
linearized compressible (example IV); and compressible (example V) flow.) 

In figure 14 lines of constant Q and 0 are plotted m the 

cp\Jr-plane. The flow direction 0 varies along the mean streamline 
(i|r = 0.5) indicating that the channel is curved. The solution is 
unsymmetrical. As for examples I and II; the lines of constant Q and 
0 are orthogonal (appendix E) . The ratio of the spacing of lines of 
constant Q and 0 is given by equation (F5) m appendix F 

( ss )q /§e\. ( 4jt/l8o\ 2* . 

(BS Te (bqJ q ( 0.03 J ^ ~ 2.7 Q 

In figure 15 lines of constant cp and ^ are plotted m the 

physical xy-plane. The shape of the channel walls is that required to 
result m the prescribed velocity distribution given by equations (35) 
and (36) and plotted m figures 2 and 13. The upstream channel width 
is twice the downstream width in order that the upstream velocity be 
half the downstream velocity. It is interesting to note that; before 
curving in the direction of the elbow turning angle; the inner wall 
first curves m the opposite direction. This behavior of the inner wall 
geometry is necessary in order to maintain the prescribed constant 
velocity along the outer wall where the velocity would otherwise 
decelerate because of the necessary curvature m the direction of elbow 
turning. This feature of the elbow geometry will also be noted in 
examples IV and V. As usual the streamlines and velocity-potential lines 
are orthogonal (appendix E) ; and; for equal increments of cp and 
form a square network (equation (F7); appendix F; with p equal to 1.0). 

In figure 16 lines of constant Q and 0 are plotted m the 
physical xy-plane. The lines of constant Q and 0 are orthogonal 
(appendix E) ; and the ratio of the spacing of lines of constant Q and 
0 is the same as that given for the same lines of constant Q and 0 
in the cfty-plane (equation (41) ) . 


24 


NACA TN 2593 


Example IV 

The fourth numerical example is the design of an elbow with the 
same prescribed velocity Q, as a function of arc length, used m 
example m but for linearized compressible flow (y = -1.0) . 

Prescribed velocity distribution . - The prescribed velocity distri- 
bution Q is the same as that for example III and with q^ equal 
to 0.80176. The variation in Q with s along one channel wall is 
plotted in figure 2. The values of q a and q-^ in equations (14a) and 

(14b) are equal to q^ and q^, or 0.40088 and 0.80176, respectively. 
For these values of qg and q^ and for the prescribed velocity 
distribution with linearized compressible flow, the elbow turning angle 
given by equation (24) was 104.08° compared with a value of 104.07° 
obtained from the relaxation solution and a value of 89.36° obtained 
for incompressible flow (example III) . 

Results . - The results of example IV are presented in figures 17 
to 19 and in tables III and IV. 

In figure 17 lines of constant q and 0 are plotted in the 
transformed cp* \|r* -plane . The solution is unsymmetncal. The lines of 
constant q and 0 are orthogonal (appendix E) , and the ratio of the 
spacing of lines of constant q and 0 is given by equation (F9) m 
appendix F. 

-f bQ \ 9 /4,t/l80\ q ft q 

X&STJ \&q/ p* - ^ 0.02 } p* ~ 0.9 p* 

where p* is related to q by equations (13) and (13b) . 

In figure 18 lines of constant cp*/A^* and ^ * are plotted m 
the physical xy-plane (where the constant A>|f* is given by equation (32) 
and is equal to 0.73782 for q^ equal to 0.80176). The shape of the 

channel walls is that required to result m the prescribed velocity 
distribution used m example III but with linearized compressible flow 
and for q^ equal to 0.80176. From continuity considerations the 
upstream channel width is 1.5385 times the downstream width. As m 
example III the inner wall of the elbow first turns in the opposite 
direction to the elbow turning angle. As usual the streamlines and 
velocity-potential lines are orthogonal (appendix E) . The ratio of the 
spacing of the lines of constant (p*/A^* and \|r*/Ailr* is given by equa- 
tion (F10) in appendix F 

j 
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(& s ) cp * (&V*\ 11/6 1 1 

p* ' 1/6 p* ' p* 

Thus the ratio of the spacing of lines of constant cp*/A\|r* and \|r*/A>|(* 
m figure 18 is a measure of the density p*. 

In figure 19 lines of constant q and 6 are plotted m the 
physical xy-plane. The lines of constant q and Q are not m general 
orthogonal (appendix E) . 


Example V 


The fifth numerical example is the design of an elbow with the same 
prescribed velocity Q, as a function of arc length, used in examples III 
and IV but for compressible flow (y = 1.4) . 

Prescribed velocity distribution . - The prescribed velocity distri- 
bution Q is the same as that for examples III and IV but with q^ 
equal to 0.79927. The variation m Q with s along one channel wall 
is plotted m figure 2. 

Results - The results of example V are presented m figures 20 and 
21 and m tables V and VI. 


In figure 20 lines of constant cp/A\|f and t|//At|t are plotted m 
the physical xy-plane (where the constant A^ is given by equation (26) 
and is equal to 0.71054 for q^ equal to 0.79927). The shape of the 
channel walls is that required to result m the prescribed velocity 
distribution used m examples III and IV but with compressible flow 
(r = 1.4) and for q^ equal to 0.79927 The upstream channel width 
is 1.5412 times the downstream width, and the turning angle is 105.31° 
compared wirh 104.07° for linearized compressible flow (example IV) and 
89.36° for incompressible flow (example III) The streamlines and 
velocity-potential lines are orthogonal, and from equation (P7) of 
appendix F the ratio of the spacing of the lines of constant cp/A^r and 
iJr/Aty is given by 


(&S)tp 

T&ST7 


1 _ 1/6 1 
p “ 1/6 p 


by 

&cp, 


1 

P 
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Thus, as for linearized compressible flow (example IV), the ratio of the 
spacing of lines of constant cd/a^t and i|r/A\|f m figure 20 for equal 
increments of Cp/Avjf and is a measure of the density p. 

The shape of the elbow for compressible flow (example V, fig. 20) 
is nearly the same as the shape of the elbow for linearized compressible 
flow (example IV, fig. 18). Therefore, m figure 21 the contours of the 
walls for both examples are compared. The difference m contours is 
Very small and it is concluded that, if a nonviscous gas with arbitrary 
y (1.4, for example) were to flow through a channel designed for linear- 
ized compressible flow (y = -1.0), the resulting velocity distribution 
along the channel walls would be nearly the velocity distribution pre- 
scribed for the linearized compressible flow, at least if the linearized 
flow were selected (by the choice of q a and q-^) so that the densities 
were equal for both types of flow at the maximum and minimum velocities 
and if the ratio of these prescribed velocities is not too large (2:1 m 
the numerical examples) . This conclusion is important because the design 
method for linearized compressible flow is considerably faster than the ' 
design method for compressible flow with y other than -1.0. 


SUMMARY OF RESULTS AND CONCLUSIONS 

A general method of design is developed for two-dimensional 
unbranched channels with prescribed velocities as a function of arc 
length along the channel walls. The method is developed for both com- 
pressible and incompressible, irrotational, nonviscous flow and applies 
to the design of elbows, diffusers, nozzles, and so forth. Two types 
of compressible flow are considered: the general type with arbitrary 

value for the ratio of specific heats y (1.4, for example) and the 
linearized type m which y is equal to -1.0. In this report (part I) 
solutions are obtained by relaxation methods on a transformed plane the 
coordinates of which are the streamlines and velocity-potential lines 
m the physical plane; m part II solutions are obtained by a Green 1 s 
function. The present method of solution gives complete information 
concerning the flow throughout the channel. 

Five numerical examples are given ana the results presented m 
plots of lines of constant velocity and flow direction or lines of 
constant physical coordinates m the transformed plane and streamlines 
and velocity-potential lines or lines of constant velocity and flow 
direction in the physical plane. Among the five examples are three 
elbow designs for the same prescribed velocity as a function of arc 
length along the channel walls but with incompressible, linearized 
compressible, and compressible flow. The numerical results of these 
three elbow designs are tabulated to enable a detailed comparison of 
the three designs. 
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The shapes of the elbows for compressible flow and for linearized 
compressible flow are very nearly the same and it is concluded that, if 
a nonviscous gas with arbitrary y (1.4, for example) were to flow 
through a channel designed for linearized compressible flow (y = -1.0), 
the resulting velocity distribution along the channel walls would be 
nearly the velocity distribution prescribed for the linearized com- 
pressible flow, at least if the linearized flow were selected so that 
the densities are equal for both types of flow at the maximum and 
minimum velocities and if the ratio of these velocities is not too large 
(2:1 m the numerical examples). This conclusion is important because 
the design method for linearized compressible flow is considerably 
faster than that for compressible flow. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, July 25, 1951 
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APPENDIX A 


Symbols 


The following symbols are used in this report: 


A, B; C ; D 


A,B 


P 

Q 

q 

q : 

s 

s 


u 

x,y 


r 

6 


coefficients; equation (29) 
arbitrary constants; equation (Cla) 
coefficient; equation (14a) 

coefficient; equation (14b) 

distance m xy-plane measured normal to direction of flow 
(expressed as ratio of characteristic length equal to 
channel width downstream at infinity) 

static pressure (expressed as ratio of stagnation density 
multiplied by stagnation speed of sound squared) 

velocity .(expressed as ratio of characteristic velocity 
equal to constant channel velocity downstream at infinity) 

velocity (expressed as ratio of stagnation speed of sound) 

velocity used m linearized compressible flow and related 
to q by equation (13b) 

spacing between lines m xy- or cp\|r-plane 

distance m xy-plane measured along direction of flow (ex- 
pressed as ratio of characteristic length equal to channel 
width downstream at infinity) 

velocity parameter related to q* by equation (18) 

Cartesian coordinates m physical plane (expressed as ratios 
of characteristic length equal to channel width downstream 
at infinity) 

ratio of specific heats 

increment of 

flow direction m physical xy-plane (measured in counter- 
clockwise direction from positive x-axis) . 


0 
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A 0 channel turning angle, equation (12) 

p density (expressed as ratio of stagnation density) 

p* density used in linearized compressible flow and related to 

p by equation (13a) 

velocity potential and stream function, respectively, used 
as Cartesian coordinates m transformed plane, defined by 
equations (3) and (4) 

boundary value of \|r along left channel wall when faced m 
the direction of flow, equation (26) 

velocity potential and stream function, respectively, for 
linearized compressible flow and used as Cartesian 
coordinates m the transformed cp*i)r*-plane, defined by 
equations (15) and (16) 

A\|f* boundary value of ty*, for linearized compressible flow, 

along left channel wall when faced m the direction of 
flow, equation (32) 

Subscripts: 

ci,b quantities related to two velocities (q a and q-^, respec- 

tively) for which density given by equation (8a) is equal 
to density p given by equations (13), (13a), and (13b) 

d o conditions downstream at infinity 

Q,q along lines of constant Q and q, respectively 

u conditions upstream at infinity 

x,y along lines of constant x and y, respectively 

A \(r* left channel wall, when faced m direction of flow, along 

which \|r* is equal to A\jr* 

0 along lines of constant 6 

co,ij/,c D*,^* along lines of constant cp, \|r, r o* , and \|f*, respectively 

0 right channel wall, when faced m direction of flow, along 

which or is equal to 0 

left channel wall, when faced m direction of flow, along 
which ^ is equal to 1.0 



1.0 



30 


NACA TCT 2593 


APPENDIX B 


EQUATIONS OF CONTINUITY AND IKROTATI ONAL FLUID MOTION IN TERMS 
OF TRANSFORMED cp,^ COORDINATES 

Consider the two-dimensional irrotational motion of a fluid particle 
in the physical xy-plane. The fluid particle is defined by adjacent 
streamlines (constant \|r) and velocity-potential lines (constant co) 
spaced 5n and 5s apart as indicated m figure 22. The velocity Q 
is parallel to the streamlines and normal to the velocity-potential lines. 

Continuity . - From continuity considerations of the fluid particle 
m figure 22 


d 

3i 


(pQ 5n) = 0 


or 

d lo &e P d lQ fee Q 1_ d(5n) 

3s 5s 5n —$s 

But, from geometrical considerations (reference 6, p. 167, for example) 



_1_ d(&n) _ 30 
5n os *" 3n 

© 

and 

_i_ 5(&s) 3e 

5s dn 3s 

so that equation (Bl) becomes 

S log e p d log e Q se 

4 * ^ = 0 

os os on 

or 

S log e p 3 log e Q ^co <ii|r 

ds + 5cp ds + (5\jr dn = ® 


(B2a) 


0 


(B2b) 


which, combined with equations (3) and (4) , becomes 
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1 (d l°g e P a lQ g e Q \ be 

p \ + 5q5 ; + ^ " ° (5 ^ 

Equation (5) is the continuity equation expressed in terms of 
co,\|f coordinates • 

Ir rotational fluid motion . - For irrotational motion of the fluid 
particle in figure 22 


3H ^ Q5s) “ 0 


or 

3 loe e g + 1_ 3(Ss) 
§n 8s dn 

But, from equations (B2h) and (B3) 

s lo Se Q _ be = 

5n 5s 


(B3) 


B log e Q ^ ^0 dqp 

5\jr dri " 5q> ds, “ 0 

which, combined, with equations (3) and (4) , becomes 

a lQ g e § be „ 
p 5* 5cp = 0 


( 6 ) 


Equation (6) is the equation for irrotational fluid motion expressed m 
terms of the cp,^ coordinates. Equations (5) and (6) vere originally 
derived m modified forms by Chaplygin (reference 7) and are given m 
reference 6, page 169. 
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APPENDIX C 


RELATION BETWEEN VELOCITY AND DENSITY ASSUMING LINEAR VARIATION 
IN PRESSURE WITH SPECIFIC VOLUME 

The approximate, linear relation between pressure p and specific 
volume l/p first suggested by Chaplygin (reference 7) is given by 


from which 

dp _ _B_ 

p 2 

where A and B are arbitrary constants. 

If p denotes the static pressure expressed as a ratio of the 
stagnation density multiplied by the stagnation speed of sound squared, 
Bernoulli's equation is 


(Cla) 


(Clb) 


+ q dq = 0 
P 

which combined with equation (Clb) integrates to give the approximate 
relation between velocity and density 


-A;- - %- = constant 

2p 2 2 


For convenience equation (C2) can be written as 


1 




1 


or 


p* = (1 + 


a*V 1/2 


(13) 


where 


P* = k-jP 


(13a) 
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and 


q* = k 2 q 


(13b) 


The constants kp and kg replace the two arbitrary constants 
m equation (C2) and their values are determined so that for any two 
arbitrary values of q (designated by q a and q^) the values of p 
given by equation (13) equal the values of p given by equation (8a) . 
Thus the values of p given by equation (13) for q equal to q a or 
q^ are correct; for all other values of q the values of p 
approximate ♦ 
ditions 


a 
are 


The constants kp and kg are determined from the con- 


P* = ^pp a 


A 


q a “ k 2 q a 


Pb ” k l^b 


^b = k 2 q b 


(C3) 


J 


From equation (13) and the conditions given by equation (03) 





1 - 


1 - 


'a a 

Pb% 


% 


(14a) 


and 


ko = — 



(14b) 


where p a and are determined by equation (8a) for the selected 

values of q d and q^, respectively. 
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The values of ^ and q ^ mighty for example, he selected to equal 
the maximum and minimum values of q (which values of q must occur on 
the channel walls and are therefore known) . Also, the values of q a and 

q^ might he selected to equal the upstream and downstream velocities 
q^ and q^. In this case the upstream and downstream channel widths 

would then satisfy continuity for a gas with the correct value of y 
(1.4, for example) . If the upstream and downstream velocities are equal, 
their value and the value of some other velocity (the maximum or minimum 
velocity, for example) can he selected for q a and q^ or, if desired, 
q a can he equal to q-^ m which case if 

q a = q + € where e — > 0 

q b = <1 

it can he shown from equations (14a) and (l4h) that 

(C4a) 

and 




This latter case, m which q a = % = <l> corresponds to the method 
used by Chaplygin (reference 7) and Karman-Tsien (reference 9) m which 


the correct relation between 


and — is replaced by a straight line 


(equation (Cla)) that is tangent to the correct relation at one point 
(where = %) . 
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APPENDIX D 


EQUATIONS OF CONTINUITY AND IRROTATIONAL FLUID MOTION IN 
TERMS OF TRANSFORMED cpf,i|r* COORDINATES 

Consider the two -dimensional irrotational motion of a fluid particle 
m the physical xy-plane. The fluid particle is defined by adjacent 
streamlines (constant ^*) an< ^ velocity- potential lines (constant cp*) 
spaced 6n and 6s apart as indicated m figure 22. The velocity q* 
is parallel to the streamlines and normal to the velocity-potential lines. 

Continuity . - From continuity considerations of the fluid particle 
m figure 22 


(p*q*Sn) = o 


or 


5 lQ g e P* d lQg e q* -L 5 (Sn) 

5s + 5s + Sn os - ® 

which combined with equation (B2a) becomes 


p " dcp* d lo & e q* dcp * dt* 


d log^ 

dcp* ds + dcp* ds + d^* dn ~ ® 


or, from equations (15) and (16) 


■L /5 log e p* 5 logg q*\ 


"5cp 




dcp ! 


5^ 


= o 


But, from equation (13) 


q * 2 5 log e q* 


! 5 log e p* 
^ d<? * * 


so that equation (Dl) becomes 


l 


1 d lp ge q* 50 = 0 

^2 dcp* + 5^* 


1 + q 


(Dl) 


(D2) 
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Finally, if 


q* 

1 +/ 1 + q * 2 


(18) 


then 


d log e q* 

/ ? = & l0 S e u 

Vi + i * 2 

bo that equation (D2) becomes 

B log e u 

— + W" 0 


(D3) 


(17) 


Equation (17) is the continuity equation expressed in terms of 
coordinates and log e 


Irrotational fluid motion . - For irrotational motion of the fluid 
particle in figure 22 




s - ° 


or 


8 log e g 1 8(6s) 
5n + 6s on 

which combined with equation (B2b) becomes 


0 


d lQ ge dj|r* be dcp* _ 

dijr * dn dcp* ds 

or, from equations (13), (15), and (16) 

1 d log e q* be _ n 
I p cty* dtp* 

VI + q* 

Finally, from equations (D3) and (D4) 


(D4) 


S log u 
e 

ch|r* 



( 20 ) 


Equation (20) is the equation for irrotational fluid motion expressed in 
terms of cp*,i|r* coordinates and log e u. 
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APPENDIX E 


ORTHOGONAL CURVES IN q*|r- AND xy- PLANES 

If, for example, lines of constant Q and 9 are orthogonal m 
the cp\)r-plane the product of their tangents equals -1.0. This condition 
is satisfied if 


SQ dQ d0 _ n 

Sep Sp SJjr 


But, from equations (5) and (6) 


(El) 


Sq be Sq be _ 

Sep dep S\jf S\jF 


\(1 f 

\ de be 

s log e p be - 

L(p 

j 3q> " 

dep cK0_ 


(E2) 


so that for compressible flow equation (El) is not, in general, satisfied 
and therefore lines of constant Q (or q) and 9 are not orthogonal 
in the cpijr-plane. For incompressible flow p is equal to 1.0 and the 
right side of equation (E2) is zero so that equation (El) is satisfied 
and therefore lines of constant Q and 0 are orthogonal m the 
cpijr-plane . 


From equations (lla) to (lid) in differential form 


Bx _ cos 0 _ By 
cfcp~ Q ~ p chjr 


By _ sin 0 _ Bx 
cKp ” Q ~ ^ 


so that 


bx by bx by _ /, 2s Sx by , _s 

+ (E3j 

For compressible flow the right side of equation (E3) is not, m general, 
zero so that lines of constant x and y are not orthogonal m the 
cp\(r-plane. For incompressible flow the right side of equation (E3) 
becomes zero so that lines of constant x and y are orthogonal m the 
cp^r-plane . 
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From the usual definitions of cp and \|r 


d<P 

5x 


. - 1 5\|r 

= Q cos 0 = — ■£-*- 
P oy 


dcp a -1 5i|r 

= Q sin 0 = — ^ 




so that 


Sep cty , Sep S\|r ’ 

+ = 0 

Thus, for both compressible and incompressible flow lines of constant 
cp and aJt are orthogonal in the xy-plane. 

In terms of Q and 0 the equations for continuity and lrrota- 
tional motion in the xy-plane reduce to 


s log e Q 

Sx 


S0 


sm 0 cos 0 


S log e p 

5y 


2 a iog e P 

cos 0 ^ 


S log e Q 


"Sy" 


S0 

SE 


sin 0 cos 0 


S log e p 

Sx 


- sm 0 


S log e p 

5y 


so that 



Sq S0 _ 


Sx 5x 3y Sy 

„ 50 

f 5 iog e p 

' Q SE 

(sin 0 cos 9 ^ + 

50 

/ d log e P 


Ism 0 cos 0 ^ + i 


S log e 

dx 

d log e p 

Sy - 


(E4) 


For compressible flow the right side of equation (E4) is not, m general, 
zero so that lines of constant Q (or q) and 0 are not orthogonal 
m the xy-plane. For incompressible flow the right side of equation (E4) 
becomes zero so that lines of constant Q and 0 are orthogonal m the 
cpnjr-plane. 
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Likewise for linearized compressible flow it can be shown that 


dq dO dq d0 _ 
dcp* dcp* + d\jr * dcp* “ 0 


(E5) 


dx by bx by 

dcp* dq>* dijr* dijr* 


= (1 - 


r> *2'v dx by 

P ’ dip chp 


(E6) 


and 


dq>* dp ^ 

^T5T + 


da>* by* _ 
° 


dq d0 dq d0 

S S ^y <5y 


(E7) 


d0 i Q Q 

q ^ |sm 0 cos 0 


d log e p* 


e h 2 * lo Se P* 

+ cos 6 ^ 


- q ^ (sin 0 cos 0 


d log e p* 

SE 


+ sin^ 0 


d log e p* 

~5y 


) 


(E8) 


Thus, from equation (E5) lines of constant q and 0 are orthogonal 
m the cp*ty*-plane and from equation (E7) lines of constant cp* and 
p are orthogonal in the xy-plane. But from equation (E6) lines of 
constant x and y are not orthogonal m the cp*p-plane and from 
equation (E8) lines of constant q and 0 are not orthogonal m the 
xy-plane . 
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APPENDIX F 


RATIO OF CURVE SPACING FOR SETS OF ORTHOGONAL CURVES 
IN cpty- AND xy- PLANES 

Consider for example the case of orthogonal lines of constant Q 
and 0 in the cp\|r-plane (incompressible flow, appendix E) . If (dS) g 

is the differential distance along a line of constant 0 between two 
curves of constant Q 


(dS)/ = (dCp)/ + (d+j/ (FI) 

where the subscripts 0 indicate that the changes are made along a line 
of constant 0. The change m Q along (dS)g is 


dQ = ^0 + 5 * ( F2 ^ 

Also, because &6 is zero along (dS)^ 

0 - 3 $ We + If (a,,) 0 ( p3 > 


From equations (Fl) to (F3) 


(dS) 




Likewise, if (dS)^ is the differential distance along a line of 
constant Q between two curves of constant 0 


(F4a) 


(<3 S) 


2 

Q 


(d0) 2 


(HT - (^T 

dQ S0\ 2 
\dcp di|r ” d\|r Sep/ 


(F4b) 
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Thus, from equations (F4a) and (F4b) the ratio of curve spacing for 
orthogonal lines of constant Q and 0 m the cp^-plane becomes 



which, from equations (5) and (6) with p equal to 1.0, becomes 



Likewise it can be shown that for incompressible flow m the 
Cfnjr-plane with lines of constant x and y 


(5S) X ^ ^ 


(F6) 


For both compressible and incompressible flow in the xy-plane with lines 
of constant cp and \|f 

/6£\ i 

For incompressible flow m the xy-plane with lines of constant Q and 0 


(F7) 


(8S) q 

IST7 



(F8) 


For linearized compressible flow m the cp*ty*-plane with lines of 
constant q and 0 


(6S) q (* 9 \ q 

WT e \&q / p* 


(F9) 


And for linearized compressible flow m the xy-plane with lines of 
constant cp* and 


( &S V / S\Jr * \ 1 

(8s) t * " [wj p* 


(F10) 
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TABLE I - DISTRIBUTION OF VELOCITY Q AND FLOW DIRECTION 0 Hi TRANSFORMED cp\J/ -PLANE FOR EXAMPLE III (ELBOW WITH INCOMPRESSIBLE FLOW) 
[Prescribe! variation In Q with arc length s along channel walls plotted in fig 2 , ^ = 05, = 1 0, A0 = 09 36° ] 
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TABLE III - DISTRIBUTION OF VELOCITY q AND FLOW DIRECTION B IN TRANSFOFMED -PLANE 

FOR EXAMPLE IV (ELBOW WITH LINEARIZED COMPRESSIBLE FLOW) 


[Prescribed, variation in Q with arc length 8 along channel walls plotted in fig 2, 

0=05, Q, =10, q = 0 80176. = 0 73782, A 6 = 104 07° ] 

u. a a 



0 

1/6 

1/3 

1/2 

2/3 

5/6 

1 

0 



q 

• 

q 

6 

q 

0 

q 

6 

q 

• 

q 

• 

q 

9 

- 11/6 

0 4009 

0 


0 4009 

0 


0 4009 

0 


0 4009 

0 


0 4009 

0 


0 4009 

0 


0 4009 

0 


- io/s 

4009 


01 

4009 


01 

4009 


00 

4009 


00 

4009 


00 

4009 

- 

01 

4009 

- 

01 

- 9/6 

4009 


01 

4009 


01 

4010 


01 

4010 


00 

4010 

- 

01 

4009 

- 

01 

4009 

- 

01 

- 8/6 

4009 


02 

4010 


02 

4010 


01 

4010 


00 

4010 

- 

01 

4010 

- 

02 

4009 

- 

02 

- 7/6 

4009 


03 

.4010 


03 

4011 


01 

4011 


00 

4011 

- 

01 

.4010 

- 

.02 

4009 

- 

03 

- 6/6 

4009 


05 

4011 


04 

4012 


02 

4013 


00 

4012 

- 

02 

4011 

- 

04 

4009 

- 

05 

- 5/6 

4009 


08 

4012 


07 

4015 


04 

4015 


00 

4014 

- 

04 

4012 

- 

07 

4009 

- 

08 

- 4/6 

4009 


14 

4015 


12 

4019 


07 

4020 


00 

401 B 

- 

07 

4014 

- 

11 

4009 

- 

13 

- 3/6 

4009 


24 

4019 


20 

4025 


11 

4027 

- 

01 

4024 

- 

11 

4017 

- 

18 

4009 

- 

21 

- 2/6 

4009 


40 

4026 


33 

4037 


1 ? 

4039 

- 

03 

4033 

- 

20 

.4022 

- 

31 

4009 

- 

35 

- 1/6 

4009 


70 

4041 


57 

4057 


27 

405 C 

- 

05 

4048 

- 

33 

4030 

- 

50 

4009 

- 

56 

0 

4009 

1 

31 

4070 


96 

4093 


36 

4090 

- 

17 

4071 

- 

58 

4042 

- 

84 

4009 

- 

92 

1/6 

4072 

2 

82 

4141 

1 

49 

4155 


43 

4139 

- 

40 

4104 

- 

99 

4058 

-1 

34 

4009 

-1 

45 

2/6 

4243 

3 

88 

4268 

1 

77 

4251 


24 

4207 

- 

87 

4148 

-1 

63 

4080 

-2 

07 

4009 

-2 

22 

3/6 

4489 

4 

00 

4444 

1 

46 

4377 

- 

38 

4295 

-1 

70 

4202 

-2 

59 

4106 

-3 

11 

4009 

-3 

28 

4/6 

4780 

3 

17 

4654 


50 

4526 

-1 

49 

4396 

-2 

92 

4265 

-3 

90 

4135 

-4 

46 

4009 

-4 

65 

5/6 

5094 

1 

44 

4882 

-1 

17 

4689 

-3 

16 

4506 

-4 

62 

4333 

-5 

63 

4167 

-6 

21 

4009 

-6 

41 

6/6 

5415 

_ 

98 

5118 

-3 

43 

4857 

-5 

34 

4621 

-6 

76 

4403 

-7 

75 

4200 

-8 

33 

• 4009 

-8 

52 

7/6 

5732 

-4 

00 

5352 

-6 

23 

5024 

-8 

01 

4734 

-9 

35 

4472 

-10 

29 

4232 

-10 

85 

4009 

-11 

04 

8/6 

6036 

-7 

48 

5578 

-9 

48 

5186 

-11 

10 

484 4 

-12 

34 

4539 

-13 

22 

4262 

-13 

74 

4009 

-13 

91 

9/6 

632 C 

-11 

35 

5793 

-13 

12 

5340 

-14 

57 

4947 

-15 

70 

4601 

-lb 

49 

4291 

-16 

97 

4009 

-17 

13 

10/6 

6602 

-15 

52 

5994 

-17 

09 

5482 

-18 

37 

5043 

-19 

30 

4659 

-20 

09 

4317 

-20 

52 

4009 

-20 

67 

11/6 

6855 

-19 

96 

61/8 

-21 

33 

5613 

-22 

46 

5131 

-23 

.34 

4712 

-23 

98 

4341 

-24 

37 

4009 

-24 

49 

12/6 

7086 

-24 

62 

6346 

-25 

80 

5732 

-26 

79 

5210 

-27 

56 

4759 

-28 

12 

4362 

-28 

46 

4009 

-28 

58 

13/6 

.7293 

-29 

46 

6496 

-30 

47 

5837 

-31 

32 

5281 

-32 

00 

4801 

-32 

49 

4381 

-32 

79 

4009 

-32 

89 

14 /6 

7477 

-34 

45 

.6629 

-35 

31 

5931 

-36 

03 

5343 

-36 

62 

4838 

-37 

05 

4398 

-37 

31 

4009 

-37 

40 

15/6 

7636 

-39 

56 

6744 

-40 

27 

6013 

-40 

89 

5398 

-41 

39 

4872 

-41 

77 

4413 

-42 

00 

4009 

-42 

08 

16/6 

7769 

-44 

/6 

6842 

-45 

34 

6084 

-45 

86 

5447 

-46 

30 

4902 

-46 

64 

4427 

-46 

65 

4009 

-46 

93 

17/6 

7875 

-50 

02 

.6924 

-50 

47 

6146 

-50 

90 

5492 

-51 

30 

4931 

-51 

63 

4440 

-51 

85 

.4009 

-51 

93 

18/6 

7953 

-55 

27 

6991 

-55 

61 

6202 

-55 

98 

5537 

-56 

36 

4961 

-56 

72 

.4456 

-56 

98 

4009 

-57 

08 

19/6 

8001 

-60 

.49 

7045 

-60 

72 

6257 

-61 

06 

.5586 

-61 

47 

4999 

-61 

91 

4477 

-62 

28 

4009 

-62 

44 

20/6 

8018 

-65 

55 

7091 

-65 

68 

6315 

-66 

03 

5647 

-66 

53 

5055 

-67 

13 

.4515 

-67 

79 

4009 

-68 

16 

21/6 

8018 

-70 

32 

7135 

-70 

45 

6385 

-70 

85 

5730 

-71 

49 

5142 

-72 

37 

4598 

-73 

46 

4072 

-74 

80 

* 2/6 

8018 

-74 

81 

7186 

-74 

96 

6471 

-75 

43 

5840 1 

-76 

22 

5272 

-77 

36 

4745 

-78 

91 

4243 

-81 

03 

23/6 

8018 

-78 

97 

7245 

-79 

15 

6573 

-79 

68 

5978 

-80 

59 

5441 

-81 

93 

.4948 

-83 

79 

4489 

-06 

33 

24/6 

8018 

- 82 , 

.82 

7311 

-83 

01 

5690 

-83 

59 

. 6138 ; 

-84 

57 

5641 

-86 

02 

5190 

-88 

01 

4780 

-90 

69 

25/6 

8018 

-86 

28 

7381 

- 86 , 

.48 

6816 

-87 

07 

.6313 

-88 

08 

5860 

-89 

55 

5455 

-91 

55 

5094 

-94 

16 

26/6 

8018 

-89 

37 

7452 

-89 

56 

6947 

-90 

15 

6494 

-91 

14 

6089 

-92 

57 

5729 

-94 

48 

.5415 

-96 

93 

27/6 

8018 

-92 

07 

7523 

-92 

25 

7077 

-92 

81 

6676 ; 

-93 

75 , 

631 C 

-95 

10 

6002 

-96 

88 

5732 

-99 

11 

28/6 

.8018 

- 94 . 

.40 

7591 

-94 

57 

7203 

-95 

09 

68521 

-95 

97 

6540 

- 97 , 

.21 

.6268 

-98 

83 

6038 

-100 

83 

29/6 

.8018 

-9b 

39 

7654 

-96 

54 

7321 

-97 

02 

| 7020 ! 

-97 

82 1 

6753 

-98 

94 

6522 

-100 

39 

6329 

-102 

17 

30/6 

8018 

-9e 

05 

7713 

-98 

20 

7432 

-98 

62 

7177 ! 

-99 

34 

6952 

-100 

34 

6759 

-101 

63 

6602 

-103 

19 

31/6 

8018 

-99 

44 

*766 

-99 

56 

7532 

-99 

94 

' 7321 

-100 

58 

7135 

-101 

46 

6978 

-102 

59 

6855 

-103 

95 

32/6 

8018 

- 100 , 

.56 

7813 

-100 

67 

7623 

-101 

01 

: 7451 

-101 

56 i 

7301 

-102 

33 

7178 

-103 

31 

7086 

-104 

49 

33/6 

8018 j 

-101 

47 

.7855 

-101 

57 

7702 

-101 

85 

j 7566 

-102 

33 

7449 

-102 

99 

7357 

-103 

83 

7293 

-104 

84 

34/6 

8018 

-102 

18 

7890 

-102 

26 

7772 

-102 

51 

! 7667 

-102 

91 

7580 

-103 

47 

7515 

-104 

18 

| 7477 

-105 

03 

35/6 

8018 

-102 

73 

7921 

- 102 , 

.80 

7831 

! -103 

00 

7753 

-103 

34 

7691 

-103 

80 

i 7651 

-104 

39 

7636 

-105 

10 

36/6 

8018 

- 103 , 

.14 

7946 

-103 

20 

7880 

-103 

37 

7825 . 

-103 

64 , 

7785 ; 

-104 

01 

7764 

-104 

48 

7769 

-105 

05 

37/6 

8018 ; 

-103 

45 

79661 

-103 

49 

7920 

-103 

62 

! 7883 

-103 

83 

7860 ; 

-104 

12 

7855 

-104 

49 

7875 

-104 

91 

38/6 

8018 ; 

-103 

66 

7982 

-103 

69 

7950 

-103 

79 

.7927 

-103 

95 

7917 

-104 

16 

7923 

-104 

42 

7953 

-104 

71 

39/6 

8018 

-103 

81 

7°94 

-103 

83 

7973 

-103 

90 

1 7960 

-104 

02 

7957 ; 

-104 

16 

7968 

-104 

33 

8001 

-104 

49 

40/6 

8018 

-103 

90 

8002 

-103 

92 

7989 

-103 

97 

, 7902 

-104 

05 

7982 

-104 

14 

! 7993 

-104 

23 

8018 

-104 

28 

41/6 ! 

8018 

-103 

97 

8008 

-103 

98 

8000 

-104 

01 

7996 

-104 

06 

7997 

-104 

11 

1 8004 

-104 

16 

8018 

-104 

18 

42/6 

8018 1 

-104 

00 

8011 1 

-104 

01 

8006 

-104 

03 

8004 

-104 

06 

8005 ! 

-104 

09 

8010 

-104 

12 

8018 

|-104 

13 

43/6 

8018 

-104 

03 

8014 

-104 

03 

8011 

-104 

05 

8009 

-104 

06 

8010 

-104 

08 

8013 

-104 

10 

8018 

i -104 

10 

44 /e 

8018 

-104 

04 

8015 

-104 

05 

8013 

-104 

05 

8013 

-104 

06 

8013 

-104 

08 

8015 

-104 

08 

8018 

-104 

09 

45/6 

. 8018 , 

-104 

05 

8016 

-104 

05 

.8015 

-104 

06 

8015 ; 

-104 

07 

8015 

-104 

07 

8016 

-104 

08 

8018 

-104 

08 

46/6 

8018 

-104 

06 

8017 : 

-104 

06 

8016 

-104 

06 

8016 ; 

-104 

07 

8016 

-104 

07 

8017 

-104 

07 

8018 

-104 

07 

47/6 

8018 

-104 

06 

8017 1 

- 104 , 

.06 

8017 

-104 

06 

8016 : 

-104 

07 

8017 

-104 

07 

8017 

-104 

07 

8018 

-104 

07 

48/6 

8018 

-104 

06 

6017 , 

-104 

06 1 

8017 

-104 

06 

8017 : 

-104 

07 

8017 

-104 

07 

8017 

-104 

07 

! 8010 

-104 

07 

49/6 

8018 ' 

-104 

06 

8017 

-104 

06 

8017 

-104 

06 

8017 

-104 

07 

8017 

-104 

07 

1 8017 

-104 

07 

8018 

-104 

07 

50/6 

8018 

-104 

.06 

8018 

-104 

06 

8018 

-104 

06 1 

8010 

-104 

07 

8018 

-104 

07 

8018 

- 104.07 

8018 

-104 

07 



46 


NACA TN 2593 


TABLE XV - DISTRIBUTION OF PHYSICAL COORDINATES x AKD y Hi TRANSFORMED -PLANE FOR 

EXAMPLE IV (ELBOW WITH LINEARIZED COMPRESSIBLE FLOW) 


[prescribed variation in Q with arc length s along channel walla plotted in fig 2, 

Q = 0 5, Q =10, q = 0 80176, M* = 0 73782, A 0 = 104 07° ] 
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TABLE V - DISTRIBUTION OF VELOCITY q AND FLOW DIRECTION 6 IN TRANSFORMED -PLANE FOR 
EXAMPLE V {ELBOW WITH COMPRESSIBLE FLOW {7=1 4)) 


[ Prescribed variation In Q vith arc length e along channel walls plotted in 
fig 2, = 05, = 1 0, q a = ° 79927, 0 71054, A0 = 105 31° ] 
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Figure 7 - Lines of constant velocity Q and flow direction Q in physical xy-plane for example I 

Incompressible flow, prescribed velocity given in figure 2 
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Figuie 8. 
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- Prescribed velocity distribution as function of arc length along channel wall for example II. 

Equation (38) . 


90£Z 


NACA TN 2593 






NACA TN 2593 


59 


x 



X 


Figure 11. - Streamlines and velocity-potential lines in physical xy-plane for example II. 
Incompressible flow, prescribed velocity given m figure 8. 
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Figure 12 - Lines of constant velocity Q and flow direction 0 m physical xy-plane for 

example II. Incompressible flow, prescribed velocity given in figure 8. 





Figure 13. - Prescribed distribution of log e Q as function of 9 along channel 

walls for examples III, IV, and V. 
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Figure 15. - Streamlines and velocity-potential lines m physical xy-plane for example III 
Incompressible flow; prescribed velocity given in figures 2 and 13 
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Figure 16. - Lines of constant velocity Q and flow direction 6 m physical xy-plane for 
example III. Incompressible flow, prescribed velocity given in figures 2 and 13 
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Figure 18. - Streamlines and velocity-potential lines in physical xy-plane for example IV 
Linearized compressible flow, prescribed velocity as function of arc length along channel 
walls same as for example III (fig. 2) and with equal to 0.80176. 
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Figure 19. - Lines of constant velocity q and flow direction 0 in physical xy-plane for 
example IV. Linearized congress lble flow, prescribed velocity as function of arc length 
along channel walls same as for example III (fig. 2) and with equal to 0 80176. 
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Figure 20 - Streamlines and velocity-potential lines in physical xy-plane for example V 

Compressible flow (y = 1.4), presented velocity as function of arc length along 
channel walls same as for examples III and IV (fig. 2) but with equal to 0,79927. 
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Figure 21 - Comparison of channel wall shapes for compressible flow (example V) with y equal to 1 4 

and for linearized compressible flow (example IV) for same prescribed velocity as function of arc 
length along channel walls (fig. 2) 
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